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ABSTRACT 

An experimental investigation has been carried out on the heat transfer 
characteristics of a turbulent, negatively buoyant, two-dimensional wall jet. 
This flow is encountered in many practical circumstances such as those related 
to heat rejection and enclosure fires. However, very little information exists 
in the literature on this problem. A jet of hot air is discharged adjacent to a 
water-cooled, isothermal surface in an extensive environment. The heat transfer 
to the surface from the jet, whose discharge temperature is taken as greater 
than the surface temperature, is measured as a function of distance along the 
isothermal surface for several values of wall, jet and ambient temperatures. 

The conditions at the jet inlet are expressed in terms of the Grashof number Gr 
and the Reynolds number Re, wide ranges of which are investigated. The total 
heat transfer to the isothermal Bt ecr is obtained by integrating the variation 
of the local heat flux along the surface. 

It is found that the total heat transfer rate to the isothermal. vertical 
surface decreases with an increase in the mixed convection or buoyancy parameter 
Gr/Re- The effect of the surface temperature if on the downward penetration of 
the jet and on the heat transfer rate is also investigated. The penetration 
distance is found to decrease with an increase in T.- Correlating equations are 
derived from the experimental data to characterize the dependence of the surface 
heat transfer rate on the various physical parameters that govern the transport 
process. An interesting flow circumstance arises when the surface temperature 
is higher than the ambient temperature but lower than the jet discharge 
temperature, since the downward wall jet flow interacts with an upward buoyancy- 
induced flow adjacent to the wall. The heat transfer mechanisms are 
investigated in detail for this flow. Several other interesting aspects are 


observed and investigated. 


AT 


NOMENCLATURE 


cross sectional area of the slot through which the jet is discharged 
specific heat of air at constant pressure 

width of the slot for jet discharge 

magnitude of gravitational acceleration 

Grashof number, Gr = gA(T,-T,)D°/v 

local surface heat transfer coefficient 

height of the enclosure 

thermal conductivity of air 

height of the isothermal vertical surface 

mass flow rate in the upward natural convection flow for the 
case when Ts aacle 

mass flow rate discharged by jet per unit length of the slot 
momentum flow rate per unit slot length at the jet inlet 
Nusselt number based on D, Nu, = hD/k 


Nusselt number based on 7 Nu. = ens 


local heat transfer flux 

dimensionless local heat transfer flux, q = q"D/kAT 

total net heat transfer to the isothermal surface from by the 
heated jet flow 

total thermal energy input by the jet Qin ~ Epencoron ours 
Reynolds number, Re = UD/v 

Richardson number, Ri = gA(T-T,)D/U- = Gr/Re" 

local fluid temperature 


surface temperature excess over the ambient, AT = th og aie 


T, discharge temperature of the jet 

Ts temperature of the isothermal surface 

od temperature of the surroundings 

U, discharge velocity of the jet 

ae : maximum velocity in the upward natural convection flow adjacent to 


the plate for the case of i Pat 


W width of the vertical plate 

x vertical coordinate distance, measured downward from slot 

B coefficient of thermal expansion of the fluid 

v kinematic viscosity of the fluid 

Fe penetration distance of the jet, measured downward from slot 

by vertical distance up to which heat is transferred from the jet to 


the isothermal plate 
p- density of air at the jet discharge 
6 dimensionless temperature, @ = Oh Bee VAG Ua tae, 


7] dimensionless surface temperature §@ = (T_-T_)/(T_-T_) 
s 5 2 on 


INTRODUCTION 
Buoyant jets are frequently encountered in nature as well as in industry. 
In nature, they are important in meteorological, oceanographic and environmental 
studies. The rejection of thermal energy and of chemical waste products to the 
atmosphere and to water bodies involves turbulent buoyant jets, such as the 
flows emerging from cooling towers or chimneys. Heat extraction and energy 
storage systems, such as those employed in solar energy utilization, are also 


often concerned with buoyant jet flows [1,2]. 


Laminar buoyant jets have been analyzed in several investigations. These 
studies clarified the basic mechanisms underlying such flows. For instance, 
Mollendorf and Gebhart [3] carried out a perturbation analysis for a vertical, 
ania axisymmetric jet, with a small amount of thermal buoyancy. They found 
that, as expected, an aiding buoyancy effect increases the velocity and an 
opposing one decreases it. The buoyancy effect was found to increase with 
decreasing Prandtl eee due to an increase in the ellie ndoi thickness of 
the thermal boundary layer. The laminar buoyant jet flow, with small buoyancy 
effects, was also studied by Schneider and Potsch [4]. They obtained the first 
approximation to the effect of buoyancy on the flow by the method of matched 
asymptotic expansions, for 0.5 < Pr < 1.5. Trends similar to those obtained 
from the analysis of Mollendorf and Gebhart [3] were observed and the range of 
validity of the analysis was clarified. Several other studies have obtained 
numerical and experimental results on laminar buoyant jets, as reviewed by 
Jaluria [5]. However, turbulent jets are of much greater practical interest and 


have, therefore, been studied much more extensively [1]. 


Very little work has been done on jets in which the buoyancy effects oppose 
the flow, such as heated jets discharged downward. Such flows arise in room 
fires and in mixed convection transport in enclosures. Energy extraction and 
heat rejection processes often lead to regions where the jet flow-is opposed by 


the buoyancy force, see, for instance, Refs. [6,7]. 


The negatively buoyant, axisymmetric, jet has been studied experimentally 
and analytically, using integral models, as reported by Turner [8] and Seban, 
Behnia and Abreu [9]. They found that the Flow penetrates to a fintee height, 
where the vertical velocity drops to zero and flow reversal occurs, as expected. 
They also reported that the flow spreads outward horizontally as it approaches 
the location of flow reversal. The analysis did not consider the effect of the 
reverse flow on the jet, which is itself driven by an external momentum source. 
The idealized circumstances of a point-source jet with momentum and buoyancy 
input, but zero mass flux at the inlet, was considered in the analysis. 
Penetration distances were found to be close to the experimentally determined 
values. But the flow rates and the velocity and temperature distributions were 


not accurately predicted by the analysis. 


In this paper, the heat transfer to an isothermal vertical surface in the 
mixed convection circumstance of a negatively buoyant, turbulent, wall jet 
flowing adjacent to the surface is considered in detail. As mentioned above, 
this problem is of particular interest in room fires and in energy extraction 


and storage systems. 


The problem of combined forced and natural convection transport over a 


vertical surface has been studied by several researchers [10]. Among the 


earliest were Sparrow and Gregg [11] and Szewczyk [12], who obtained a series 
solution using perturbation from the basic natural and forced convection laminar 
flows. Both aiding and opposing flows were considered and solutions were 
obtained for velocity and temperature profiles, shear stress and heat transfer 
rates, over a wide range of Pr and of the perturbation parameters. Merkin [13] 
employed a parameter expansion, rather than a coordinate expansion. For 

the circumstance of the buoyancy force opposing the laminar forced flow, a 
series solution valid near the leading edge was obtained. The results obtained 


by Merkin [13] are of particular relevance to the present work. 


There are many practical situations where turbulent, negatively buoyant, 
wall flows arise. For example, when the fire plume in an enclosure hits the 
ceiling, it spreads out and finally turns downward at the corners [14,15]. At 
this stage, the flow of the wall jet is downward, whereas the buoyancy force is 
upward, since the flow is at a temperature higher than that of the surroundings. 
Similarly, heated flows rising adjacent to two walls may flow along the ceiling, 
meet at the top and be pushed downward, against the buoyancy force. This 
results in a negatively-buoyant, two-dimensional, free jet. At a later stage in 
the growth of an enclosure fire, a hot upper layer is established above a 
relatively cooler lower layer. Then, a buoyancy-induced flow arises in the 
lower layer, adjacent to the wall which is usually hotter than the neighboromg 
fluid. This flow becomes negatively buoyant as it penetrates into the upper 
layer, across the interface separating the two regions [16,17]. Negatively 
buoyant flows may also result due to the discharge of fluid at a temperature 
different from the ambient air and water media in which the transport occurs. 
Such flows arise, for instance, in energy extraction from a solar pond and heat 


rejection with jets inclined downward [7,18]. 


A downward negatively buoyant flow is retarded due to the buoyancy force, 
which acts vertically upward. This results in a decrease in the flow velocity, 
which ultimately becomes zero and then reverses the direction. Therefore, the 
negatively buoyant flow penetrates to a finite depth in the ambient environment 
before the reversal occurs. It also entrains fluid from the ambient mediun, 
resulting in a larger flow rate at the level of the source, as compared to that 


discharged by the jet itself. 


Goldman and Jaluria [19] carried out a detailed experimental investigation 
to study the basic flow characteristics of negatively buoyant wall and free 
jets. An insulated surface was employed in the study on wall jets. The 
penetration distance $5 was measured and found to be largely governed by the 
mixed convection or buoyancy parameter Gr/Re“, which is also often termed the 
Richardson number Ri. All these dimensionless parameters are based on the inlet 
conditions of the jet, as defined later. The net entrainment into the flow was 
also determined. It was found that as Gr/Re? increases, resulting in larger 
buoyancy effects, the flow penetration decreases and a stronger reverse flow, in 
the direction of the buoyancy force, arises. The entrainment into the flow was 
found to increase with Ger eae over the range considered, as a result of this 
increased reverse flow at larger Gr/Re- 

The effect of the thermal conditions at the wall, i.e., whether it is 
adiabatic, heated or cooled, on the behavior of negatively buoyant flows has not 
been investigated in detail. Also, since the wall was adiabatic in the 
experimental work of Goldman and Jaluria [19], no data on the resulting heat 
transfer was obtained. The dependence of the wall flow on the thermal 


conditions at the surface and the heat transfer to the wall are of considerable 


importance in the modelling of the transient heating of the walls in enclosure 


fires and in the other practical cases mentioned above. 


It is clear from the above review that the heat transfer characteristics of 
a turbulent, negatively buoyant, two-dimensional, wall jet is not available in 
the literature. The objective of the present study was to investigate such 
flows and to determine the heat transfer from the flow to a cooled, isothermal, 
surface. Heat flux data have been obtained for several wall temperatures, 
varying Gr/Re*from 0.05 to 1.05. The local heat flux distribution has been 
integrated over the surface up to oo! the jet penetration distance to obtain the 
total heat transfer from the jet to the surface. The effect of the wall 
temperature on the penetration of the jet has also been studied. It is found 
that, for a given ambient temperature, the penetration distance $5 decreases as 
the surface temperature Ts increases, for Gr/Re” > 0.1. For smaller values of 
Gr/Re*, es was found to be essentially independent of the wall temperature. The 
effect of the upward natural convection boundary layer flow, which arises when 
the surface temperature = is greater than the ambient temperature T_, is also 
studied and related to the basic transport mechanisms in the flow. Several very 
interesting and important trends are observed and considered in terms of the 


flow characteristics of a negatively buoyant wall jet. 


Lr 


EXPERIMENT 

Experimental Arrangement 

Figure 1 shows a sketch of the experimental arrangement. A blower sends 
ambient air, over a fairly wide range of flow rates, through a copper tube which 
is 5 cm in diameter and 1 min length. The copper tube is heated by means of 
three fiberglass insulated strip heaters wrapped around it. The energy input to 
each of the heaters is varied by means of power controllers. A diffuser at the 
end of the copper tube helps in discharging the heated air as a two-dimensional 
jet. The width of the slot for discharge may be varied up to about 0.1m. 
Several diffuser designs were considered in order to ensure a uniform two- 
dimensional flow at the exit. The temperature and velocity measurements 
indicated that fairly uniform distributions, with variations of less than 10%, 
were obtained at the diffuser exit. The heat losses from the copper tube and 
from the diffuser to the environment were minimized by employing insulation made 


of foam and fiberglass. 


The discharge velocity of the jet Uy could be varied from about 0.3 to 2.5 
m/s. The discharge temperature ba could be raised to as high as approximately 
150°C. Thus, this arrangement could be used for studying fairly wide ranges of 
the governing parameters Re and Gr, based on the inlet conditions, as defined 
later. Values as high as around 3500 and 10° for Re and Gr, respectively, 


could be obtained [19]. 


The heated, two-dimensional, air jet is discharged into a glass enclosure, 
1.5 m high and 1.37 m x 0.3 m in cross section. The bottom of the enclosure is 
kept open in order to allow the wall jet flow to entrain ambient air from below 


or to flow out of the enclosure at a relatively small velocity, resulting in a 


BL 


negligible effect of the outflow on the negatively buoyant flows under 
investigation. The top and the far side of the enclosure are also kept open to 
allow a free exchange with the ambient medium. This arrangement, thus, 
simulates a wall jet in an extensive ambient medium and, as seen later, ina 


steady flow circumstance. 


A water cooled aluminum plate was located vertically adjacent to the side 
wall of the enclosure over which the heated jet was discharged. The plate was 
91.4 cm in height. The remaining portion of the side wall was insulated. The 
details of the plate design are shown in Fig. 2. Four rectangular copper tubes, 
each 2.5 cm x 1.2 cm in cross section, run along the length of the plate. The 
tubes were kept flush with the plate surface and were bonded to it by soldering. 
It was ascertained that the contact between the two surfaces was good by means 
of temperature measurements. As can be seen from this figure, water circulates 
over approximately 44% of the plate area. Water from an external tank enters at 
the top of the plate, with that emerging from this arrangement being allowed to 
drain into a sink. The temperature of the water entering the plate assembly 
could be maintained at a desired value by mixing hot and cold water streams from 
two separate tanks. Nine thermocouples were embedded in the plate to monitor 
the plate temperature. Details on the thermocouple assembly for plate 


temperature measurement is also shown in Fig. 2. 


Several 0.3 cm diameter blind holes were drilled into the plate. 
Thermocouples, supported by 0.3 cm diameter glass tubes, were inserted into 
these holes. It was ensured that the thermocouple beads always touched the 
plate surface and that a good contact was maintained, so that the plate 


temperature was accurately obtained. Several similar arrangements have been 


Tat 


used by the authors in the past and this background was used in the design of 
the arrangement. A 2 cm plastic cube was attached to the plate to support each 
glass tube. Finally, the plate was insulated at the back. It was found that 
the plate could be maintained at a constant temperature indefinitely. Also, it 
was confirmed by measurement that a fairly uniform temperature distribution was 
obtained at the plate surface. The maximum temperature difference measured 


between any two surface thermocouples was around 0.5°C. 


The temperature of the heated jet at the inlet was measured using a rake of 
five thermocouples. The rake was placed just below the slot through which the 
jet is discharged. The average of these five temperatures was taken as the jet 
temperature Ty in the experiments. The temperature was found to be very 
uniform, within +2°C, across the slot cross-section. A DISA constant 
temperature hot-wire anemometer was employed for the velocity measurements. The 
hot wire was calibrated by measuring the frequency of the vortices behind a long 
cylinder [20]. The frequency in the wake allows the determination of air 
velocity across the cylinder. The velocity distribution at the discharge was 
again found to be very uniform, within variation of less than 5%. The average 
of this distribution was taken to represent the jet discharge velocity U- 

Thus, fairly uniform velocity and temperature distributions were obtained at the 
discharge, and the measurements yielded the corresponding values of the jet 


discharge velocity U, and temperature ae: 


The heat transfer to the plate was measured by means of microfoil heat flow 
sensors (RdF Type 20472-3). Each of the heat flow sensors was 15 mm x 6 mm in 
cross-section and 1 mm in thickness. They could easily be attached to the plate 


surface. The electric output (in millivolts) from each heat flow sensor was 


13% 


converted into the corresponding heat flux q" (in ayn) with the help of 
individual ceushee wien curves supplied by the manufacturer. The calibration was 
also checked for accuracy, using electrically heated foils. These heat flow 
sensors were found to be highly sensitive and accurate. They covered the entire 
range of heat flux obtained in the experiments. Eighteen such heat flow sensors 
were attached, at a regular intervals of 5 cm, along the centerline of the 
plate. The sensors were also moved around to confirm the two-dimensionality of 


the transport process. 


The voltage signals from the hot wire anemometer, thermocouples and the 
heat flow sensors were fed into a data acquisition system. The signal was also 
recorded on a strip chart recorder. The fluctuating component was generally 
found to be small, being less than 5% of the mean temperature and heat flux 
values. Thus, the mean values, which are of main interest in Pirie cuays were 
obtained easily. The experimental arrangement simulates a hecactcele buoyant 
wall jet in an extensive isothermal environment. Considerable care had to be 
taken to ensure that a high level of repeatability, within 5-10%, was maintained 
and that accurate measurements of the tempreature field and of the heat transfer 


rate were obtained. 
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RESULTS AND DISCUSSION 

The width D of the long slot discharging the jet is generally taken as the 
characteristic dimension for a jet in an extensive environment [9,19]. Then, 
the dimensionless parameters that govern the flow are the Reynolds number Re, 
the Grashof number Gr and the Richardson number Ri, which is also often known as 
the mixed convection parameter, Gr/Re. The heat transfer at the wall may be 
presented in terms of the Nusselt number Nu, [10]. These dimensionless 
quantities are defined for a negatively buoyant wall jet in an extensive 


isothermal medium at temperature T, as: 


3 
U_D eBCE =TsyD 
ROM he ee 
Re arcs Cr ™ Ct 
gB(T,-T)D hD 
Ri = Sheene aye a and Nu, = eat (2) 
Oo 


where Uy and T, are the velocity and temperature at the discharge of the jet, g 
is the magnitude of the gravitational acceleration, 6 is the coefficient of 
thermal expansion of the fluid, h is the local convective heat transfer 
coefficient and v is the kinematic viscosity of air at the discharge temperature 
TS: 


It must be mentioned that the above dimensionless parameters are obtained 


when the governing equations and the boundary conditions are nondimensionalized 


with US» T, and D as the characteristic quantities [10]. Another Nusselt 
number, Nu,, based on the penetration distance on is also defined as 
hé 
Nu. = (3) 


6 k 


Fo: 


where o is the penetration depth of the thermal effects in the wall jet flow, 


as defined quantitatively in terms of the thermal field later. If the average 
heat transfer coefficient h is employed instead of h in Eqns. (2) and (3), the 


average Nusselt number Nu, or Nu. is obtained. The dimensionless local 


temperature § and the dimensionless surface temperature a are defined as 


T-TL La oe 
bm aE 5 Test (4) 
1) «a oO foe) 


where T is the local temperature in the flow and T. is the surface temperature. 
The surface heat flux q" is appropriately nondimensionalized in terms of the 


above characteristic quantities to yield the dimensionless heat flux q as 


a= 15)- 


The data consist mainly of the temperature distributions in the wall jet 
flow and of the local heat flux measurements along the isothermal surface for 
different values of the physical variables such as temperatures and flow rates. 
Detailed measurements were taken for several values of the wall temperature TS. 
The mixed convection parameter Ri = Gr/Re- was varied from about 0.05 to 1.05 in 
the experiments, by varying the discharge velocity Us slot width D and 
temperature T.° However, for convenience, the slot width D is held constant at 
6.5 cm for the results presented here. The effects of varying US» i. and D were 
all found to be well correlated in terms of the Richardson number Ri, as found 
by Goldman and Jaluria [19]. Also, even though D is used as a characteristic 
dimension here, results may be presented in terms of the discharge flow rate, 


per unit slot length, m,, momentum inflow rate and energy input rate [5,8]. For 


16. 


instance, m/pM., where Me is the momentum input per unit slot length at the 


jet discharge, may be used as a characteristic dimension instead of D. 


Figure 3 shows the vertical temperature profiles taken at y/D = 0.154 for 
the negatively buoyant jet flow at Ri = Gr/Re* = 0.123 and 0.263. This value of 
y corresponds closely to the location of maximum temperature in the wall jet. 
Thus, Fig. 3 indictes the decay of the maximum temperature in the flow with 
vertical distance away from the jet inlet, x = 0. It is seen that the 
temperature level decreases sharply along the plate up to a certain distance. 
This is followed by a more gradual decrease, with the temperature finally 
becoming equal to that of the surroundings. The penetration distance Pe 
represents the penetration of thermal effects in the flow. This penetration is 
indicated quantitatively by a sharp rise in the temperature level as one 
proceeds from the bottom to the top of the tank. It was found that this 
location was quite sharply defined, within 1-2 cm, and is marked in the figure. 
A similar behavior has been buseread earlier [19] and repeatability to within 
5-10% of the value of hs was obtained in separate experiments. The vertical 
temperature distribution shown in Fig. 3 is important in the understanding of 
the basic heat transfer and fluid flow characteristics of a negatively buoyant 


wall jet and will be referred to later. 


The variation of the penetration distance by with the mixed convection 
parameter Gr/Re“ for adiabatic and isothermal wall conditions, considering 
three different typical values of AT, the wall temperature excess over the 
ambient temperature T_, are shown in Fig. 4, with Pe normalized by D. As can be 
seen from this figure, the penetration distance + for an adiabatic wall 


condition is very close to that reported earlier by Goldman and Jaluria [19]. 


17, 


It is also observed that penetration is essentially the same for various wall 
temperatures if Gr/Re“ < 0.1. However, for Gr/Re” > 0.1, the penetration 
distance was found to vary significantly with the wall temperature excess AT. 
For a given value of Gr/Re“, the jet penetration distance $5 was found to be 
smaller for a larger value of the wall temperature excess AT. This is an 
expected trend, since.a larger wall temperature excess AT results in a smaller 
heat transfer to the wall. This gives rise to a larger buoyancy effect, which 
retards the flow more rapidly, resulting in smaller Ay: Similar trends were 


observed for other wall temperatures considered in this study. 


The effect of the wall temperature on the penetration of the jet also 
depends on whether T; is larger or smaller than T,. When the surface 
temperature ive is higher than that of the surroundings, the surface induces a 
buoyancy driven flow which originates at the bottom edge of the surface and 
proceeds vertically upward. This convective mass flow opposes the jet flow 
which is discharged at the top. The calculations for the resulting laminar 
natural convection flow [10] for a surface temperature of 48°C with T. = 25°C, 
at the locations where the jet was found to reverse direction, show that the 
maximum natural convection velocity ee is 0.087 m/s and 0.345 m/s, for 
Gr/Re” = 0.1 and 1.05, respectively. Thus, pyle varies from 0.19 at 
cr/Re* = O21" ton0asyat eryRes = 1.05. However, the ratio of the mass flow 
Myc in the natural convection boundary layer to the discharged jet mass flow mo, 
Meo /M, > varies from 0.004 at Gr/Re" wiQbl tOn0t2), ac Gripe = = 1.05. These 
calculations indicate that the natural convective flow may have a_ considerable 

ys 


effect on the reversal of the jet flow for Gr/Re > 0.1, if the plate 


temperature T, is much larger than T|. However, in many problems of practical 


S 
interest, such as the enclosure fire at the early stages of the fire, T is 


18; 


quite close to T. and the upward natural convection flow is negligible. This 
study investigates the circumstances of both negligible and significant natural 


convection boundary layer flow along the plate. 


The variation of the non-dimensional penetration distance Dh as a 
| 2 2 
function of the parameter Gr/Re over the range 0 < Gr/Re < 1.0 may be 
expressed in terms of appropriate correlations. These were derived from the 


data obtained for -5°C < AT < 30°C and are given as: 


6 


os = 4.5 |cx/Re* 


-0.402 
for the adiabatic wall condition (6) 


8p 2<-0 362 
Dar (2.14 - 8.466 ge) +( 2, OP tees o Os) (Gr/Re ) 


for the isothermal wall condition (7) 


where AT = To-T, and 0 = (To-T)/(T,-T,) - These correlations were found to 
represent the experimental data closely. The correlation coefficients were 
larger than 0.95, indicating this close agreement. In these correlations, the 
effect of the inlet conditions of the jet is largely represented by the 
parameter Gr/Re* and that of the wall temperature by Os: The penetration is 
also seen to be larger for the adiabatic condition, as compared to the 
isothermal cases, over the ranges considered. This is an expected result, 
since the buoyancy effect is the largest in this case ae nee of negligible 


energy loss to the insulated surface. 


Figure 5 shows the variation of the dimensionless local heat transfer flux 


q = q"D/kAT along the surface for Gr /Re- = 0.265 and AT = -3°C. It is seen that 


ee 


the local heat transfer rate decreases sharply along the plate up to x/D = 5.5 
and then remains essentially constant, at a fairly small value of q over the 
remaining portion of the surface. Thus, the local heat transfer to the surface 
does not become zero even after the jet reverses its direction. This is 
obviously due to the fact that the ambient temperature T, is slightly higher 


than the wall temperature T This temperature difference results in heat 


s* ; 
transfer by natural convection to the surface. The penetration distance oo? 
obtained from the vertical temperature distributions for these experimental 
conditions, has also been marked in the figure. It is interesting to note that 
the jet loses most of its thermal energy to the surface within a fairly short 
distance from the slot at these conditions. This AL: represents the region 
over which the heated jet loses thermal energy to the plate. It is denoted by 
by and is also shown in Fig. 5. Typically, bys was found to be 70-80% of the 
corresponding penetration distance as in the present experiments. Thus, even 
though the thermal effects, in terms of temperature disturbances and temperature 
excess over the ambient, penetrate to a distance oe the heat flux becomes 
negligibly small upstream of the location of flow reversal. Such trends have 


been observed in other convective flows as well, such as deviation from the 


laminar trends in transition to turbulence [10]. 


Figures 6 and 7 show the variation of the local heat transfer rate q with 
the vertical distance x/D for wall temperature excess AT values of 13°C and 
23°C, respectively. The basic trends are similar to those observed in Fig. 5, 
except that the heat transfer rate q becomes negative after a certain distance, 
implying heat transfer from the surface to the flow. This distance is bi as 
defined earlier, and indicates the region over which the surface gains energy 


from the flow. The corresponding values of by have also been marked on the 
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respective vertical temperature profiles in Fig. 3. As seen in Fig. 3, bs 
indicates the vertical location where the local jet temperature essentially 
becomes equal to the surface temperature To. Thus, from Figs. 5-7, it is seen 
that the negatively buoyant jet loses energy to the surface up to a distance by 
where the local temperature level in the wall jet becomes essentially equal to 
the surface temperature. The temperature in the flow from by to SS is actually 
lower than the surface ESiperature for positive AT. Consequently, heat is 
transferred from the isothermal surface to the jet in this region. ‘This is 
shown as Q_ in Figs. 6 and 7 and has wy, subtracted from Qy) also shown in the 


figure, to obtain the total net energy transfer to the plate Q. 


As can be seen from these figures, the local heat transfer q remains 
essentially constant beyond the penetration distance Fe to the bottom edge of 
the surface. The height of the surface extends up to x/D = 14.3 and, in the 
cases considered, the penetration distance 2p was much less than the height L of 
the surface. The heat flux was generally not measured to the bottom edge of the 
surface. Since the surface temperature was higher than that of the surroundings 
in Figs. 6 and 7, natural convection heat transfer arises from the surface to 
the surroundings. The measured heat flux value is expected to decrease and 
ultimately become zero near the bottom edge of the surface. This was confirmed 
in a few cases and the experimental heat transfer data were found to be close to 
that predicted by analysis of natural convection from a vertical isothermal 


surface [10]. 


The local heat flux distribution along the isothermal surface may also be 
expressed in terms of the Nusselt number Nu, based on the width D of the slot 


through which the jet is discharged. The distribution of the Nusselt number 


zie 


Nu, = hD/k along the isothermal surface, at Gr/Re“ = 0.47 and for the three 

plate temperatures considered earlier, is shown in Fig. 8. The corresponding 
penetration distances o, have also been marked in the figure. Obviously, the 
trends are the same as those observed in Figs. 5-7, since D and k are taken as 
constant for these experimental conditions. Other wall temperatures were also 


considered and similar trends were observed. 


The variation of the non-dimensional distance for heat loss, 6./D, with 
Gr/pat! for different wall temperatures, is shown in Fig. 9. It is seen that 5 
decreases as Gr /Re- increases, which is expected since a larger value of Gr/Re- 
implies larger buoyancy effects and thus a shorter penetration distance [19]. 
These results may be compared with the corresponding results for the penetration 


distance o shown in Fig. 4. It is seen from Fig. 9 that 5,, is different for 


H 
the different wall temperature excess values, even for Gr/Re-, fess Chara le 
This was not found to be the case for the Penge actor distance one see Fig. 4. 
From these results, it may be concluded that for small values of Gr/Re“, even 
though the penetration distance °5 is essentially unaffected by the wall 
temperature excess AT, by does vary significantly with AT. For large values of 


Gr/Re-, both the penetration distance es and the heat loss distance by are 


significantly dependent on the plate temperature excess AT. 


For various values of Gr/Re“ and different wall conditions, the local heat 
flux q" was obtained as a function of distance x along the plate. Then, the 
total net heat transfer Q to the plate by the jet was obtained by integrating q" 
up to the corresponding penetration distance S and multiplying the result by 
the plate width W. For wall temperatures larger than T,, the thermal energy 


transferred by the isothermal plate to the jet, Q_, was subtracted from the‘ 
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thermal energy transferred to the plate by the jet Q, to obtain Q as Q = Fag 


Q.. This gives the net energy transferred by the jet to the wall. 


The total net heat transfer Q to the isothermal plate by the jet for the 
three plate temperature excess values considered earlier is shown as a function 
of Gr/Re? in Fig. 10. It is seen that for a AT of -3°C, Q sen eam sharply up 
to Gr/Re” = 0.25. This is followed by a gradual decrease up to Gr/Re” = 1.05. 
For plate temperature excess of 13°C and 23°C, the variation in Q with Gr/Re” is 
found to be much smaller. In this figure, Q is nondimensionalized by the total 


energy input by the jet QIN which is calculated as follows 


Qin = Fences i ae (8) 
Typically, 3-15% of the jet input energy is transferred to the plate. At large 
Gr/Re“, VQ is found to be weakly dependent on Gr/Re“, indicating that the 
2 


penetration distance approaches an essentially constant value as Gr/Re 


increases to large values. This effect was also seen in Fig. 4. 


The corresponding results in terms of the average Nusselt number Nu, = 


hD/k, where h = SUMO and Wé, is the heat transfer area,are shown in Fig. 11. 


The trends are similar to those observed in Fig. 10. It is interesting to note 


that the average Nusselt number Nu, tends to increase to large values as Gr/Re” 
becomes small and approaches zero at large Gr/Re’. The trends are clearly 
linked with the dependence of the penetration distance on Gr/Re” . At small 
Gr/Re, on is large and the jet loses much of its energy to the surface, whereas 


at large Gr/Re“, ae is small and only a small portion of the jet energy is lost 


Lhe 


to the surface. Also, at a given Re, the inlet temperature Ts is lower at a 
smaller Gr/Re?, resulting in a smaller heat transfer rate to the surface. Thus, 
there are two competing effects. As Gr/Re“ becomes smaller, the area for heat 


transfer increases but the temperature difference (T5-TS) goes down. The 


combination of these opposing effects also gives rise to a maximum in Nu, as 


seen in Fig. 11. 


An average Nusselt number Nu, based on the penetration distance De may also 


be defined as 
Shs koe pee GG (9) 


since h = eee oT - The variation of this average Nusselt number Nu. with Gr/Re 
is shown in Fig. 12 for different wall conditions. It is again seen that the 
Nusselt number decreases sharply up to around Gr/Re” <i? deli Sie. Se tnen 


followed by a gradual decrease as Gr/Re“ increases. Since the effect of S on 


the Nusselt number is eliminated in this formulation, the variation of Nu, is 


similar to that of Q, as shown in Fig. 10. 


The average Nusselt number Nu, varies with the wall temperature excess AT 


and the mixed convection parameter Gr/Re“, as shown in Fig. 11. Thus, a 


correlating equation may be derived from the data obtained to yield Nu, as a 
function of the dimensionless surface temperature excess vise and Gr/Re*. The 


result obtained is given as 


oie 


Nu, - (9.0 - 10.66 4.) (CE/Resye 2 (10) 


This equation was found to provide a good representation of the data, over the 
experimental ranges considered. A correlation coefficient of about 0.9 was 


computed, indicating the accuracy of this correlation. 
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CONCLUSIONS 

An experimental study has been carried out to investigate the heat transfer 
characteristics of a negatively buoyant, two-dimensional, wall jet. Negatively 
buoyant flows often arise in many problems of practical interest, such as 
enclosure fires, thermal discharge into the environment and thermal energy 
storage systems. A heated, two-dimensional, jet is discharged vertically 
downward adjacent to an isothermal surface in a large enclosure. A water cooled 
aluminum plate was employed and maintained at desired wall temperatures. The 
measurements consisted of the local heat flux distribution over the isothermal 
surface for various values of the mixed convection parameter Gr/Re’, varying 
from 0.05 to 1.05. The thermal field in the flow was also measured in detail. 


The data were obtained for various wall temperatures. 


The total heat transfer to the isothermal surface was determined by 
integrating the local heat flux q up to jet penetration distance oh along the 
surface. It was found that, typically, the total heat transfer to the surface Q 
iecraased with an increase in the parameter Gr/Re. This is related to the 
decrease in jet penetration as Gr/Re? is increased. It was also found that, 
with an increase in the wall temperature excess over the ambient, the total heat 
transfer to the surface Q decreases, as expected. The effect of wall 
temperature on jet penetration has also been investigated in detail. It is 
found that the jet penetrates to a shorter distance when the wall is kept at a 
temperature higher than that of the surroundings, at a constant value of Gr/Re’. 
This is attributed to the fact that, when the wall is warmer than the ambient, 
it induces a vertical upward, natural convection boundary layer flow adjacent to 


the wall. This flow opposes the jet which is discharged downward, resulting in 


a reduced penetration distance af the jet. 
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Experimental arrangement for the study of a heated, two-dimensional, 
downward discharged wall jet in an isothermal ambient medium 
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surface 


Typical vertical profiles in the negatively buoyant wall jet flow 


Variation of the penetration distance §6_ with the mixed convection 
parameter Gr/Re for adiabatic and threé isothermal wall conditions 


Local heat transfer flux distribution along the surface for Gr/Re* = 


0.265 at a surface temperature Te which is 3°C less than the ambient 
temperature 


Local heat transfer flux distribution along the surface for Gr/Re“ 
0.470 and a surface temperature excess over the ambient AT of 13°C 


Local heat transfer flux distribution along the surface for Grype- 
0.473 and AT = 23°C 


Variation of the,local Nusselt number Nu, along the isothermal 
surface at Gr/Re = 0.470 


Variafion of the nondimensional distance for heat loss 5,,/D with 
Gr/Re 


Variation of the net energy lost to the plate by the jet Q, 
nondimensionalized by the botal energy input Qi? with the mixed 
convection parameter Gr/Re 
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Variation of the average Nusselt number Nu, with Gr/Re 
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Variation of the average Nusselt nymber based on 6_ as the 
characteristic distance with Gr/Re P 
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